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We perform extensive MD simulations of two-dimensional systems of hard disks, focusing on the 
on-collision statistical properties. We analyze the distribution functions of velocity, free flight time 
and free path length for packing fractions ranging from the fluid to the solid phase. The behaviors 
of the mean free flight time and path length between subsequent collisions are found to drastically 
change in the coexistence phase. We show that single particle dynamical properties behave anal¬ 
ogously in collisional and continuous time representations, exhibiting apparent crossovers between 
the fluid and the solid phase. We find that, both in collisional and continuous time representation, 
the mean square displacement, velocity autocorrelation functions, intermediate scattering functions 
and self part of the van Hove function (propagator), closely reproduce the same behavior exhibited 
by the corresponding quantities in granular media, colloids and supercooled liquids close to the glass 
or jamming transition. 


I. INTRODUCTION 

A system of two-dimensional (2D) hard disks is one 
of the simplest models of a classical fluid. Despite the 
apparent simplicity, the transport properties and the na¬ 
ture of the phase transitions remain areas of active in¬ 
vestigation ever since the pioneering work of Alder and 
Wainwright mm some 50 years ago. The appearance 
of slow power-law decaying tails in the velocity auto¬ 
correlation function for moderately dense systems has 
deeply changed the classical understanding and the en¬ 
suing formulation of the kinetic theory (see for instance 
Refs. |3l l^). Moreover, controversy about the nature 
of the fluid-solid transition occurring at higher densities 
has persistently involved generations of scientists, debat¬ 
ing on whether its best representation is the Kosterlitz- 
Thouless-Halperin-Nelson-Young scenario 01 , or a first 
order transition H. This debate seems to have only 
recently arrived to a conclusion thanks to the initial mas¬ 
sive use of event-chain Monte Carlo algorithm [9], and by 
later adopting local Monte Carlo algorithm and event- 
driven molecular dynamics simulations m- 

A system may be followed via an external clock rep¬ 
resented by continuous time t , or an internal clock as¬ 
sociated with the number of collisions n of each particle. 
The continuous time representation is said to be subordi¬ 
nated to the underlying collisional process HJ. Though 
the system is clearly characterized by collision events, 
a systematic analysis of this system in collisional rep¬ 
resentation, spanning from the fluid to the solid phase, 
is missing. To the best of our knowledge the only an¬ 
alytic and numerical studies in this sense concern the 
non-Poissonian nature of the number of collisions within 
a certain time interval in the dilute limit [MU- In what 
follows we study statistical and dynamical properties of 
a 2D system of hard-disks in collisional representation, 
with packing fractions ranging from the fluid to the solid 
regimes. 


The first part of our analysis focuses on the statistics 
of on-collision velocity, free flight time r and the path 
undergone by a disk between subsequent collisions, the 
free path length (|£|). Within this context we study cor¬ 
relations between the x and y components of the velocity, 
and the dependence of free flight times on the on-collision 
velocities in systems with various packing fractions span¬ 
ning from the fluid to the solid phase. We also compare 
our findings to the kinetic approach and Enskog theory, 
reporting inconsistencies even in the fluid phase. Fur¬ 
thermore, our analysis of (r) and (|£|) as a function of 
packing fractions has unveiled plateau-like regions corre¬ 
sponding to the fluid-solid coexistence phase, which to 
the best of our knowledge have not been reported. 

The second part focuses on dynamics, comparing sin¬ 
gle particle dynamical observables in continuous time (as 
a function of t) and in collisional representation (as a 
function of the number of collisions, n). We find a full 
analogy, in spite of the fact that the process is not truly 
Poissonian [13] . We show that the velocity is never a 
Markovian process in either representations and in both 
fluid and solid phases m ■ Furthermore we have ana¬ 
lyzed the mean squared displacement, intermediate scat¬ 
tering function and self-part of the van Hove function (or 
propagator) in both representations, showing remarkable 
similarities to glassy behavior. 

We consider a 2D system of size L x x L y where L x 

( L y ) corresponds to the length in the x ( y ) direction, 

consisting of N hard disks with diameter a. The packing 

fraction r) of the system is defined as the ratio between 

2 

the area occupied by the disks, Ntt and the available 

AT 2 

area L x x L y , yielding 77 = L 13 L In collisional rep¬ 
resentation we follow a single particle counting the colli¬ 
sions it undergoes with a collision index n. For a collision 
n we now define the particle’s position x,„ , the external 
(continuous) time t n , and the on-collision velocity v n . 

We run event-driven molecular dynamics (MD) [16] 
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numerical simulations. Simulations start from a ran¬ 
dom configuration, according to the desired y, obtained 
from NVT Monte Carlo simulations. The temperature 
is kept constant by scaling the magnitude of velocities 
of each hard disk such that the kinetic energy of the 
system agrees with the equipartition theorem. The sys¬ 
tem is rectangular, where the ratio between L x and L y 
is kept at \/3. Periodic boundary conditions are applied 
in both directions and the system length is adjusted to 
keep constant the number of particles N, yielding the 
desired rj. We set fcgT = 1, a = 1. and the lengths 
L x and L y are measured in units of a. Throughout 
this work we consider the following packing fractions: 
?7 = 0.3,0.56,0.695,0.713,0.72, plotted respectively in 
black, blue, red, green and magenta. We thus span the 
different phases: fluid phase for low rj, the coexistence 
phase for 0.69 < i] < 0.723 p2I US], and above that the 
solid phase (the maximum possible packing fraction is 
given by ry = Tt/Vvl « 0.907 [19]). 

II. COLLISIONAL REPRESENTATION: 

STATISTICS 

The analysis of collisional statistics is based on three 
fundamental observables: (i) Free flight times r n = 
t n .(_i — t n , the times between consecutive collisions. The 
index n represents the number of collisions suffered by 
a particle in the system, (ii) On-collision velocity v n , 
the x and y components of the particle velocity between 
collisions, (iii) Free path length, namely the vector con¬ 
necting the particle’s position on the n-th and (n + 1)- 
th collisions: = x„ +1 — x„. We first focus on the 

statistics of these measurable observables by means of 
event-driven MD numerical simulations, identifying the 
stationary distributions and the corresponding average 
values of on-collision velocities, free flight times and free 
path lengths. 


A. Instantaneous velocity distribution function 

The distribution of instantaneous velocities of the disks 
in the continuous time representation is given by the 
Maxwell-Boltzmann equilibrium distribution function: 

4 ’ MBM= 2ihf e ~ £iT ’ (1) 

where T is the system temperature and ks is the 
Boltzmann constant. In collisional representation, how¬ 
ever, it is known that the stationary distribution of on- 
collision velocities cb rn ji(v) takes a different functional 
form HU mUl]: 
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where Iq and I\ represent the cylindrical modified Bessel 
functions of 0-th and 1-st order respectively. As shown in 
Fig. [lja) the collisional distribution exhibits wider tails, 
i.e. a higher probability for higher velocities, as compared 
to its continuous time counterpart. This discrepancy can 
be understood intuitively with the following argument: in 
continuous time representation, the velocity distribution 
is calculated across snapshots of the system at different 
intervals, therefore higher velocities will be counted less 
often since their respective free flight times are shorter 
on average, spanning across less consecutive snapshots. 
In contrast, in collisional representation each on-collision 
velocity is counted only once, upon collision. 

The expression in Eq. |2| highlights two important 
features of the velocity col isional process. The first 
is that the collisional trajectory is an isotropic pro¬ 
cess, since the velocity distribution depends only on 
its speed. The speed distribution f co u (|v|) can indeed 
be derived from Eq.([2]), yielding f coU (|v|) = 27 t0 co h(v) = 


V2k B T 


( M 2 

\2k B T 



fcoii (|v|) is plotted in Fig. Hub) against the numerical 
data: this extends to d = 2 the analytic expression of the 
speed distribution in collisional representation furnished 
for d = 3 and d = 5 in Ref. [14]. 

The second property highlighted by Eq. § is that 
the velocity components are correlated within the colli¬ 
sional process, i.e. 0 co «(v) ^ <Pcoll(v x )ip C oil(v y ) although 
Vcoii{v x ) = Vcoii{v v ) (where tp C oii(v) is the single com¬ 
ponent velocity distribution). The single component ve¬ 
locity distribution function can be obtained by integrat¬ 
ing Eq.Q over the other component, i.e. ip C oii(v x ) = 
dv y (f> co u (v). The numerical evaluation of the inte¬ 
gral expression of ip C oii(v x ) is shown in Fig.[lja), perfectly 
reproducing the numerical data. We note that Maxwell 
molecules, interacting through a pairwise potential such 
as V(|r|) ~ |r| 2 [20], exhibit identical behavior in both 
continuous time and collisional representations: they are 
both isotropic, and the x and y components are in fact 
uncorrelated, yielding the same Gaussian form of the ve¬ 
locity distribution functions as in Eq. 0 m- 

Noting that the on-collision velocity is isotropic, and 
that the x and y components are statistically identical 
(even though they are correlated), we can henceforth 
limit our analysis to the single component x (or y without 
loss of generality), i.e. v n = v x and = x n+ i — x n . 


B. Free-flight time distribution function 

The probability distribution of free-flight times P(r) 
was first analyzed in Ref.s [l3j[22j[23], where it was shown 
that it cannot be purely exponential, not even in the low 
density limit |22] . This is rooted in the fact that the char¬ 
acteristic free-flight time (or its inverse the collision fre¬ 
quency) of particles with a given velocity, (r(v)), depend 
on the speed iv| mm- To explain the non-exponential 
shape of P(t) analytically, we introduce the conditional 
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FIG. 1: (Color online) Comparison of the velocity distribu¬ 
tion in continuous time and collisional representation, (a) 
Maxwell-Boltzman distribution function of single component 
velocity for the continuous time representation, <Pmb(v), plot¬ 
ted with a blue solid line, and the collisional < p co ii{v) in a 
red solid line. Black dashed lines represent, respectively, the 
usual Maxwell-Boltzmann expression, and the theoretical es¬ 
timate obtained by numerical integration of Eq.0 over the 
y-component. (b) Speed distribution function in collisional 
representation f co u (|v|): numerics are shown as a solid red 
line, while the theoretical expression f co n (|v|) is reported in 
dashed black line. ksT =? 1, a = 1, L x and L y are measured 
in units of a. 

probability of having a free-flight time r given that the 
particle has a velocity v: P(r|v). One can then derive 
the free-flight time distribution, following mm- 

/ oo 

dvP (t|v) ( f > coii ( y ). (3) 

-oo 

The conditional probability P (r|v) is related to the vis¬ 
cosity of the system [ 22 ], and is therefore interesting to 
study. By introducing the single component conditional 
probability p(r\v) (shown in Fig. [2ja)) we can also write: 

/»00 

P(r) = / dvp(T\v)ip co u{v). (4) 
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FIG. 2: (Color online) Dependence of free flight times on ve¬ 
locity. (a) The probability distribution p(r\v) for different p 
with v = 0.0±0.05 (arbitrary units). The exponential expres¬ 
sion e (dashed lines) successfully describes the numer¬ 

ical outcomes only for the lower value of p, i.e. p = 0.3,0.56. 
(b) Average free flight time as a function of v, (r(u)), for 
several p. The dashed lines are obtained from the numerical 
evaluation of Eq. 0. 


For low 77 , it can be shown that P (t|v) = ( T (v)) e ^ 
mm, a consequence of the molecular chaos assump¬ 
tion, or strosszahlansatz , which holds in the dilute limit. 
In this limit the collisional dynamics of a tagged parti¬ 
cle is exactly a Markov process in the space of veloci¬ 
ties, where the hopping probability from velocity v' to 
velocity v, given by FF(v|v'), can be calculated analyti¬ 
cally [IS]. We find that also the single component dis¬ 
tribution function p(r\v) exhibits an exponential shape 
p(r\v) = , r K' e~ T/ ('’Wb in the dilute limit p —> 0 , as 
can be appreciated in Fig. 0a). However, though both 
the complete and single component conditional proba¬ 
bilities P(t|v) and p(r\v) are exponential in the dilute 
limit, it must be stressed that the characteristic free-flight 
time given a single component velocity, (t(u)) shown in 
Fig. 0b), differs from (r (v)). Deviation of P (r|v) from 
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a pure exponential is expected at large p due to the loss 
of the molecular chaos ansatz, and to the ensuing hop¬ 
ping dynamics which is no longer Markovian in velocity 
space. This deviation is clearly observed also in the sin¬ 
gle component p(r\v), as displayed in Fig. 0a). In any 
case substituting Eq. ([2]) in Eq. ([3 ) means that P(r) can¬ 
not be an exponential either, as shown in Fig. 0a). The 
same argument can be specified to p (t|u) and ip co ii(v ) in 
Eq. Q] . 

The characteristic free-flight time is defined as (r) = 
J 0 °° dr P{t). Substituting Eq. (0) or (0 yields (r) = 
fZo dv ( T ( v ))0coh(v) or (r) = f_ x dv (r {v))ip coU {v) ac¬ 
cordingly. (r) is displayed in Fig. B a ) for different p, 
spanning from the fluid to the solid regime. Interestingly, 
we observe a plateau associated to the coexistence phase 
consistent to a scenario where spatial regions with local¬ 
ized particles (crystallites) can be observed in between 
regions of mobile particles [7[ [HI HH [Ml [Ml IMH57] . 

The relation between the characteristic free flight time 
(r) and (t(v)) follows )13]: 

(t (v))0 co n(v) = (r)to(v). (5) 

Changing to polar coordinates and integrating over the 
angle yields M- 

(T(v))/ con (|v|) = (t)/mb (M), (6) 

We find that the same equality is valid for the single 
components, as detailed in Appendix [A| 

(t (v))ip co u(v) = (t)ipmb(v). (7) 

Substituting in Eq.(0) the integral expression of ^coii(^), 
and the numerical value of (r) (Fig0a)) leads to a semi- 
analytic estimate of (t (v )), shown in Fig. 0b) to be in 
excellent agreement with simulation data. 




FIG. 3: (Color online) Free flight time (a) and free path (b) 
probability distribution functions for several values of p. The 
exponential expression e ^ > (dashed lines in panel (a)) do 
not capture the numerical outcomes in the low as well as in 
the high p regime. 


C. Free-path length distribution function 


The free-path length r is defined as |£| = |v| r, and the 
single component obeys |£| = |t;| r. In this sub-section we 
analyze the single component free-path length distribu¬ 
tion function P(|£|), displayed in Fig. 0b) for several 
values of p. In step with the analysis made for the free- 
flight distribution function in the previous sub-section, 
we can define a conditional free-path length distribution 
function as EM- 


P (r |v) 



\ P(r |v) 

) M 


and, for the single component: 


( 8 ) 


in the low limit of p, P(|£| |v) and p(|£| |f) are both 
exponentials [13] . At the same time, a marked devia¬ 
tion from the exponential form is expected for high p 
(as shown in Fig. 0a)). The conditional mean free-path, 
given a velocity v, is given by (|£(v)|) = |v|(t(v)), and 
(l£(' y )l) = M(tCu)) for what concerns the single compo¬ 
nent reported in Fig. 0b). Expressing (|£(u)|) in terms of 
(t(u)) allows the semi-analytical evaluation of the mean 
free-path single component, in excellent agreement with 
the numerics. 

The free-path length distribution function is defined 
as: 

/ OO 

dv P (r |v) ^ co n(v), (10) 

-OO 


PMI») = /*«([|J-r)^. (9) 

Using Eqs. 0 and § together with the analysis re¬ 
ported in the previous sub-section, it can be shown that 


and has been the subject of an extensive numerical anal¬ 
ysis in Ref. [38]. In this early study it was shown that: ( i ) 
P (|£|) is nearly exponential in the fluid and solid limit; 
(ii) the kinetic theory predicts a very accurate mean free- 
path at all packing fractions, namely (|£|)fct = 557 ^; (***) 














5 






FIG. 4: (Color online) Average free flight time (a) and average 
free path (b) shown as a function of r/. It is seen that they 
both exhibit a plateau in correspondence of the coexistence 
phase 0.69 < r/ < 0.716 |17 | (blurred regions), while they 
decay monotonically in the fluid and solid phases. Insets: 
blue curves show the kinetic theory predictions (Eqs.(Bl I and 
(B31 respectively), while red curves refer to the Enskog theory 
IC^pif and p5>. 


FIG. 5: (Color online) Dependence of free path length on ve¬ 
locity . (a) The probability distribution p(|£||u) for different 
?7 for v = 0.0 ± 0.05 (arbitrary units). The exponential ex¬ 
pression -— S (|g(„)|) ( ^ (dashed lines) describes the numerical 
outcomes only for the lower value of rj, i.e. r/ = 0.3,0.56. 
(b) (|£|(t))} as a function of the velocity for several PFs. The 
dashed lines are obtained by the relation (|£(v)|) = |i’|(r(u)). 


the dimensionless free path distribution as a function of 
the rescaled free-path length, (|£|)/ct-P V exhibits 


a nearly universal behavior, independent of 77 ; ( iv ) the 
rescaled free-path length distribution does not agree well 

with the zero density scaled distribution (|£|) 0 -P f J, 

where (|£|)o is the mean free-path achieved from the En¬ 
skog theory [52] ■ In step with Eq.(10), the single compo¬ 


nent path length distribution function is defined as: 

/»oo 

P(\t\) = / dvp (|£| | v)<p coU (v). (11) 


Eqs. (|T0|) and (11) yield non-exponential free path dis¬ 


tribution functions for all 77 , once the expressions for 
</>coZi(v) and <p C oii{v) are substituted (see Fig. [3](b)). In 
the very dilute case this result coincides with the find¬ 
ing in Ref. [IS]. In Fig. §b) we plot the average free- 
path with i] ranging from the fluid to the solid phase: 


(l£l) = JTdKI |£| P(|£|). As for the average free-flight 
time, the plateau can be considered as the signature of 
the coexistence phase HZ!- In Appendix[B]the single com¬ 
ponent mean free-path furnished by the kinetic theory is 


shown to be (|£|)/ci = 27 % (Eq.(B3l) and in the inset of 


Fig. gb), one can see that this estimate reproduces just 
qualitatively the numerical mean free-path (|£|), in con¬ 
tradiction to what has been observed in Ref. [35]. As a 
consequence, no universal behavior is expected from the 
rescaled single compone nt mean free path distribution 
( see Fig 18 a)). Moreover, this rescal¬ 
ing does not provide a correct collapse of the numerical 
curves, neither when the kinetic theory estimate (|£|)fct is 
replaced by the true (|£|), as displayed in Fig. [l 8 |(b). In 
Appendix B] we also provide the mean free-path length 
esti mat e ( £|)o which arises from the Enskog theory (see 
Eq.( |B5[ )). Plotting the rescaled free-path probability dis¬ 
tribution by (|£|)o corroborates the finding of Ref. [58] : 
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the agreement of the data with the zero-density limit is 
indeed unsatisfactory (Fig. [l8|(c)). In the insets of both 
panels of Fig. [4] we report the Enskog theory prediction 
for the mean free flight time and the mean free path, 
(r)o and (|||)o resp ectively. The expressions furnished in 
Eqs. ( |B5| ) and ( |B6| ) do not describe the observed behavior 
also in the limit of very dilute systems. 

Finally, defining the mean free path as (|f|) = 

dv (|£(v)|)0 co m(v) yields, in view of Eq.® and 
Eq-0, the following relation: 

(III) = (|v|)co„t(r), (12) 


or 


(III) = (M)cont(T), 


(13) 


for the single component. The same relation in Eq. (12) 


is imposed within the framework of the kinetic theory to 
define the mean free-path (see Eq. (B2) and Ref. 1301). 
Moreover it has been rigorously validated also in one¬ 
dimensional elastic rod systems, also known as Jepsen’s 
gases m- Indeed, as in ID the velocity distribution func¬ 
tion is invariant upon elastic collisions [42], the Eq. (13) 


has been shown to hold for any type of distributions 
<p(v), not only for the Maxwell-Boltzmann distribution 
Vmb ( v ). 


III. COLLISIONAL REPRESENTATION: 
SINGLE PARTICLE DYNAMICS 




Our analysis now turns to the dynamical properties of 
the single particle, or tracer. In particular we reinterpret 
the stochastic motion of a tagged particle in positional 
space x n and velocity space v n as a global collisional pro¬ 
cess, where the continuous time t is now replaced by the 
collisional index n. Correspondingly, our analysis will 
extend to the study of the properties of r n . 

In what follows we will first consider the velocity auto¬ 
correlation function, discuss the mean squared displace¬ 
ment (MSD) and the connection among them through 
the Green-Kubo relation. We then continue to investi¬ 
gate autocorrelation functions of free flight path and free 
flight time. Lastly we complete the picture by discussing 
the intermediate scattering function and the self-part of 
the van Hove function, or the propagator. 

A. Velocity autocorrelation function 

We wish to follow the argument of Markovianity and 
its intrinsic connection to the ansatz of molecular chaos 
in dilute systems. This issue is particularly important 
and compelling since it constitutes the fundamental hy¬ 
pothesis from which the on-collision velocity distribution 
function in Eq. (|2j) is derived, and which leads to the ex¬ 
ponential shape of the free-flight time conditional proba¬ 
bility distribution P(r|v) in the low density limit. Thus, 


FIG. 6: (Color online) Normalized velocity ACF. (a) Con¬ 
tinuous time. The velocity ACF exhibits persistent memory 
effects in the low PF limit and in the high PF region. The 
predicted 1/t behavior is evident in the inset, while the an- 
tipersistent tails characterize the caging effect for high r/ (b) 
Collisional representation. The velocity appears to be non- 
Markovian also in the collisional representation, showing a 
surprisingly analogy with the continuous time counterpart. 
In the inset the 1/n behavior appears to define the low r/ sys¬ 
tems, while the antipersistent memory effects dominate the 
high rj regime. 


let us briefly recall what is the theoretical framework in 
Refs. [1 % OS] |T5] that hinges on this assumption. The 
diluteness of the system guarantees the molecular chaos 
ansatz, meaning that the evolution of the tracer velocity 
distribution function <^>(v) is governed by a linear Boltz¬ 
mann equation. This, in turn, leads to the Maxwell- 
Boltzmann distribution ^(v) in the stationary regime. 
The analysis of the Boltzmann equation provided in 
Ref. [T5] clarifies that the tracer undergoes a sequence of 
collisions, performing a Markov process in velocity space, 
whose transition rates W(v|v') are analytically calcu¬ 
lated from the linear Boltzmann equation. Hence, the 
conditional characteristic free flight time (r (v)) is ob¬ 
tained through the relation (r(v)) -1 = J dv'fF(v|v'). 
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Finally, substituting this in Eq. ([5]) yields the full ex¬ 
pression of (pcoiiiy) hr Eq. ([2]). It is important to note 
that Eqs. and § are generally valid, since they must 
hold for any 77 , not only in the dilute limit 77 —>• 0. This 

means that the ratio l ' T , v )' > is a universal function of 

d 

the speed fi(|v|), independent of the packing fraction of 
the system. In view of these simple considerations, it 
is critical to verify whether the molecular chaos and/or 
Markovian assumptions are necessary and verified con¬ 
ditions in the dilute limit. However, since we restrict 
our analysis to the dynamical properties of single par¬ 
ticles, we will focus on testing the Markovian assump¬ 
tion by studying the behavior of the velocity autocorrela¬ 
tion function (ACF) in collisional representation: (voV n ). 
Furthermore, the Markov property leads to the exponen¬ 
tial shape of P(r|v) in the low i] limit, as is reported in 
Refs. [ISOS] and shown in Fig.[2ja). 

It is instructive to compare the dependence of the ve¬ 
locity ACF (voV n ) on the collision index n to the de¬ 
pendence of its continuous time counterpart (u (0 )v(t)) 
on t. In continuous time, velocity is known to be a non- 
Markovian process for any value of 77 . The pioneering nu¬ 
merical work of Alder and Wainwright mm in the dilute 
limit has shown that (u(0 )v(t)) ~ t _1 for long times. The 
discovery of these persistent memory effects greatly in¬ 
fluenced the further development of nonequilibrium sta¬ 
tistical physics of liquid states, and it is often referred to 
as the ”2D long-time-tail problem” [3]. Several theories 
based on mode-coupling theory [33], nonequilibrium sta¬ 
tistical mechanics [44] and a kinetic approach [32 have 
been devised to explain this surprising numerical find¬ 
ing, marking, in the latter case, the birth of the “Mod¬ 
ern Era” of kinetic theory (see also the partial review 
in [3]). However, the exact form of the long-time tail is 
not fully understood j46j: while an earlier result [47] es¬ 
timated the first correction to the t~ l tail as ln(f/t 0 ), 
mode coupling theories [481 [33] and the renormaliza¬ 
tion group approach [50] have led to the prediction that 
the t" 1 decay corresponds to intermediate times, while 
at longer times the asymptotic decay takes the form of 
(t-s/lnt) 1 . Extensive numerical simulations carried out 
using a direct approach [51] , a cellular automaton lattice 
gas model [52H53] and an event-driven direct hard-disk 
simulation scheme m provided the evidence in favor of 
the (f-s/lnt )” 1 long time behavior. 

In Fig. pfa) we report the single component velocity 
ACF (u(0)u(f)) for different 77 . In the inset one can ap¬ 
preciate the asymptotic f _1 behavior for the lower 77 , 
while the accuracy of our event-driven simulations can¬ 
not clearly distinguish the predicted (t-\/hit ) -1 for larger 
times. For higher 77 , well inside the coexistence region 
and the solid phase, it is possible to see how the velocity 
ACF exhibits negative antipersistent tails. This antiper¬ 
sistence is due to the backscattering resulting from the 
“caging effect”, where particles are confined to transient 
cages, resulting in temporal anticorrelations in particle 
displacements. The cage effect has been observed in col¬ 
loidal systems close to the glass transition [55H32] . su¬ 


percooled liquids [53HT3] or granular systems close to the 
jamming transition I741l76j . and mainly connected to a 
reduction of the displacements of the particles trapped in 
the cages HO El]. Here we claim that velocity anticorre¬ 
lations might reveal the presence of solid-like structures, 
characterizing the solid-fluid coexistence phase in hard¬ 
disks systems. We note that the focus here is not on the 
continuous time behavior summarized in Fig.[b](a), which 
is meant only to provide a confirmation of the overall 
behavior in continuous time, but rather on the behav¬ 
ior of the velocity ACF in collisional representation. In 
Fig. [ 6 ](b) we report (vov n ) for several 77. The inset shows 
that in the dilute limit (vov n ) ~ n~ l for large times. 
The main panel displays the progressive approach from 
the persistent to the antipersistent regime as rj increases. 
These numerical results highlight the complete analogy 
between the asymptotic behavior exhibited by the ve¬ 
locity ACF in the continuous time and the collisional 
representations for any 77 . Most importantly, it demon¬ 
strates that the Markovian assumption for velocity in 
collisional representation is never fulfilled. In particular, 
surprisingly, it does not hold in the limit 77 —>■ 0 , revealing 
that the Markov property is not a necessary condition for 
Eq. ([ 2 ]) to be valid. Nonetheless this is puzzling, since the 
exact expression provided by the Eq. ([ 2 ]) is derived by as¬ 
suming the single particle collisional velocity process to 
be Markovian in the dilute limit. 

Finally let us compare the results obtained in 2D 
with the properties of Jepsen gas, mentioned earlier. 
In ID, the velocity process in collisional representation 
(i.e. (vov n )) is Markovian only in systems with bi- 
modal (or dichotomic) velocity distributions, i.e. <f(v) = 
\ [<5(r) — c) + 5(v + c)] where S(x) stands for the Dirac’s 
delta function. The on-collision velocity process is non- 
Markovian for systems with any other distribution func¬ 
tion, including (Pmb{v) [3T]. On the other hand, in con¬ 
tinuous time, the behavior of (u(0 )v(t)) in ID consider¬ 
ably differs from the collisional representation, in con¬ 
trast to what we have shown in 2D. For bimodal ve¬ 
locity distribution functions the continuous time veloc¬ 
ity ACF is strictly exponential [3TJ 021E5], while weak 
memory effects of the form (v(0)v(t)} ~ — 1~ 3 result 
from velocity distribution functions such as Maxwell- 
Boltzmann ‘EH si EH IM], uniform ip(v) = ^0(v + 
c)9(v — c) (75] or ip(v) = y(c 2 + 7 ' 2 ) 3 / 2 |32] . where 9{x ) 
and c represent the Heaviside’s function and a constant 
respectively. For velocity distribution functions like the 
three-modal tp{v) = fiS(v) + [5(u — c) + S(v + c)] [U, 

[75] or four-modal ip{v) = ^ [S(v — c\) + 6(v + c\)] + 
[<5(tj — c 2 ) + 5(v + c 2 )] [TO (0 < /i < 1 , ci and c 2 
constants), the memory effects are stronger: (v(0)v(t)) ~ 
—f -3 / 2 . Based on the numerical evidence, it was argued 
in Ref. 133] that the non-exponential behavior observed 
in the velocity ACF in continuous time are ascribable 
only to the non-Markovian nature of the underlying col¬ 
lisional mechanism. This conjecture seems to be fulfilled 
also in 2D. 




FIG. 7: (Color online) (a) MSD in continuous time for sev¬ 
eral r ]. The solid lines refer to the Einstein-Helfand formula 
(141 while the dashed ones are obtained by the Green-Kubo 
expression (151. They both coincide yielding, after an initial 
ballistic regime, normal diffusion for low PF and an inter¬ 
mediate slowdown (plateau) for PF in the coexistence phase 
and beyond, (b) Time-dependent transport coefficient de¬ 
fined as D(t) = f g ds(v(0)v(s)). For the lowest analyzed r/, 
rj = 0.3, D(t ) increases as ~ Vint consistent with the pre¬ 
dicted (v(0)v(t)) ~ (t\/hif) -1 [48Ft50l . For the higher PF 
shown, D(t) displays a “subdiffusive” decay in correspondence 
of the caging regime in panel (a) (see Ref. 1751 1. 


B. Mean squared displacement and free path 
autocorrelation function 

We analyze the single particle mean squared displace¬ 
ment (MSD) in collisional representation and compare 
it to its counterpart in continuous time. For the sake of 
clarity, let us start by discussing the MSD in the common 
continuous time representation. The MSD is defined as 

(Sx 2 (t)) = ([x(t)~x( 0)] 2 >. (14) 

We will be referring to Eq. ( fl4| as the Einstein-Helfand 
expression, in clear connection to the well-established 
Einstein-Helfand formula for the diffusion coefficient D = 


FIG. 8: (Color online) (a) MSD in collisional representation 
for several ?y. The solid lines refer to the Einstein-Hellfand 
formula (161 while the dashed ones are obtained by the Green- 
Kubo expression (181. They both coincide, displaying normal 
diffusion for low PF and an intermediate slowdown (plateau) 
for PF in the coexistence phase and beyond, in analogy with 
the results shown in FigFFIa). (b) Transport coefficient de¬ 
fined as D n = D n exhibits a persistent in¬ 

creasing behavior consistent with ~ \/ n In n for ?y = 0.3. This 
is in agreement to the continuous time behavior displayed 
in Figjr](b) and with the power-law tails in Fig§ For the 
higher PF, D n shows the same “subdiffusive” decay of D(t ) 
(Fig[7|b)) in correspondence of the caging regime in panel (b) 
(see the detailed analysis in appendix [Cjl. 


\ Hindoo HiH. In Fig. Ra) we report the 

outcome of our numerical simulations for several rj. Af¬ 
ter a transient ballistic regime, it is clearly seen that 
systems characterized by a lower ry exhibit normal dif¬ 
fusion behavior, whereas, for higher ry, the MSD reaches 
a plateau which eventually turns into a linear regime. 
This plateau can be considered as further evidence of 
the caging undergone by a tracer during the coexistence 
phase, in analogy with colloids [58H601 [85] . supercooled 
liquids [SHUTOUTS! 155] . and granular systems (7T1I7!)] . Due 
to the stationarity of the velocity process, the MSD can 
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be also written as 


(5x 2 (t)) = 2 f ds(v(0)v(s))(t — s), (15) 

Jo 

referred to as the Green-Kubo MSD expression, in rela¬ 
tion to the definition of the diffusion coefficient according 
to the Green-Kubo relation, D = dt(v(0)v(t)} [86] . 

Fig .[7](a) shows the excellent agreement between the MSD 
calculated according to the Einstein-Helfand expression 
from Eq. ( fT~T| ) (solid lines) and that resulting from the 
Green-Kubo formula in Eq. (151 (dashed lines). The log¬ 
arithmic tail of the velocity ACF shall lead to the weak 
divergence of the diffusion coefficient D in the long-time 
limit ( D ~ [ln(oo) + const J 1 / 2 ) [28- 25]. However, as 
stressed in Ref. m, in spite of the enormous numerical 
effort devoted to the identification of the (v(0)v(t)) log¬ 
arithmic asymptotic regime, very little work has been 
dedicated to the study of its influence on the actual value 
of the diffusion coefficient [87] [88]. Indeed, the numeri¬ 
cal simulations shown in Fig. 0 a ) seem to point out the 
absence of any correction to the linear behavior exhib¬ 
ited by the MSD, in agreement with the results reported 
in Refs. (87, 89] . However, plotting the time-dependent 
transport coefficient D(t) = \ = fg ds(v(0)v(s)), 

as shown in Fig. [7](b), shows a small but persistent in¬ 
crease within the timescale corresponding to the linear 
region exhibited by the MSD, for which one would have 
expected a con stant value. For rj = 0.3 we observe 
D(t) ~ -y/ln (t), consistent with the velocity ACF form 

/ _X -1 

(v(O)v(t)) ~ I t\/\nt 1 , and in contradiction to the nu¬ 

merical finding already present in one of the earliest nu¬ 
merical works on this subject [881. For larger 77 , D(t ) ex¬ 
hibits a “subdiffusive” decay consistent with the appear¬ 
ance of negative tails of (v(0)v(t)) (shown in Fig. [ 6 ](a)) 
and with the corresponding plateau observed in the MSD 
(shown in Fig. [7][a)) for intermediate times; asymptot¬ 
ically it approaches a constant value, in line with the 
asymptotic diffusive regime. The rigorous connection be¬ 
tween anomalous diffusion and power-law tails exhibited 
by ('c(O)u(f)) is provided by the generalization of Kubo’s 
theorem in Ref. m- 

In collisional representation the analog of the Einstein- 
Hellfand MSD can be expressed as 


(Sxl) = ([z« -a:o] ) 


(16) 


and is reported in Figj8|a) (solid lines). For low 77 the 
MSD undergoes normal diffusion in the collision index 
n. We note that in collisional representation the initial 
ballistic regime is suppressed by definition. Increasing 77, 
the same intermediate plateau as in continuous time is 
observed. Since the single-component trajectory is ex¬ 
pressed as 


n— 1 


- x o = 53 


(17) 


1=0 



FIG. 9: (Color online) Normalized free path length ACF. The 
non-Markovianity of the free path length is apparent both 
at low and high 77 . For low 77 (inset) the 1/n behaviour is 
clearly shown, consistent with the velocity ACF in Fig. [ 6 ](b) . 
The systems with high 77 are displayed in the main panel, 
exhibiting a negative part, i.e. antipersistence due to the 
caging phenomenon. 


and assuming the stationarity of the free-path process, 
one arrives at the Green-Kubo expression in collisional 
representation: 


(5xl) = n(t 2 ) + 2 53 “ m )> 


(18) 


where (£ 2 ) = d££ 2 P(£). Hence, the numerical eval¬ 

uation of Eq. ( p~8] > requires the calculation of the free- 
path ACF (£o£n73 which is reported in Fig. [9j The ob¬ 
served behavior resembles that of the velocity ACF. In¬ 
deed, for any 77 the free-path process £„ is non-Markovian, 
with persistent asymptotic tails oc 71 -1 for small 77 (see 
the inset of Fig. [ 9 ]), and antipersistent as 77 increases. 
When calculated according to the Green-Kubo expres¬ 
sion in Eq. ( [T8| , the MSD (dashed lines in Fig. |s](a)) is 
shown to coincide with the Einstein-Hellfand expression 
in Eq. (161. In analogy with the continuous time analysis, 
the asymptotic nT 1 regime should lead to the divergence 


of the MSD in Eq. ([18]). However, no appreciable devia¬ 
tion from the linear MSD trend is detected, albeit in the 
range for which (£o£n) displays the n _1 behavior. We 


now recast the Green-Kubo MSD expression in Eq. (18) 


by introducing the “symmetrized” free path autocorrela¬ 
tion function as C^(n) = (£o£™) - \ (£, 2 ) 5 n ,0 [781 ISO] : 


n— 1 

i Sx l) = 2 53 C^(m)(n-m). (19) 

m —0 

We now identify the on-collision transport coefficient 
D n = X)m=oQf( m )i an d report it in FigjsVb) . This 
quantity dictates the 71 -dependence of the M5D in colli¬ 
sional representation as D(t) does for (5x 2 (t)), detailed in 
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FIG. 10: (Color online) Normalized free flight time ACF, 
C T t(ti) = (totVi) — (r) 2 . The non-Markovianity of the free 
flight time process r n characterizes both low and high 17 . 


the theorem in Appendix [Cj This theorem can be viewed 
as the extension to impact dynamics of the generalized 
Kubo theorem valid in continuous time m- In particu¬ 
lar when D n attains a constant value, one expects nor¬ 
mal diffusion (&r 2 ) ~ ra as displayed in Figj 8 ](a) for low 77 . 
However, D n behaves as ~ y/n\nn for the lowest packing 
fraction 77 = 0.3, in complete analogy with the continuous 
time analysis reported in FigjTjb), and in partial agree¬ 
ment with the observed persistent tails exhibited by the 
free path ACF (£ 0 £„) (Fig|9j). When D n exhibits a decay 
in the form of a power-law, ( Sx 2 ) is sublinear in the col¬ 
lision index n. This is indeed the case for intermediate n 
with large 77 , ultimately ascribable to the caging effect. 


C. Free flight time autocorrelation function and 
free flight and free path cross correlation function 


Since in collisional representation space and time can 
be studied separately, we now focus on the free flight time 
autocorrelation function defined as C TT (n) = (tqTu) — 
(r) 2 , reported in Fig. 10 for different 77 . As for the 


free path £ n , the free flight time r n is a stationary pro¬ 
cess. For an uncorrelated process one would expect 
(ro r n ) — (t ) 2 = [(t 2 ) — (r) 2 ] 5 n fi. However, the analysis 
in Sec. |II B demonstrates that the time elapsed between 
subsequent collisions is velocity-dependent; this velocity, 
on the other hand, is a correlated process both at low 
and high packing fractions, as discussed in Sec. me 
Therefore correlations should also characterize the dis¬ 
crete dynamics of free flight times of single particle m, 
as indeed confirmed by the numerical outcomes in Fig. |10[ 
The expression for the time ACF C TT (n) is formally de¬ 


fined through the correlation 


(t 0 t„) = 


dr 0 ■ ■ ■ dr n dv 0 • • • dv n t 0 t„ 


p(ro, ■ ■ ■ T n \v 0 , ■ ■ ■ v n )<j) co ii(v 0 , • • • v n ). 


( 20 ) 


Hence as long as times depend on velocities, the uncorre¬ 
lation of the velocity process implies the uncorrelation of 
the free flight times collisional dynamics and vice versa. 
Indeed, if we assume the on-collision velocities to be un¬ 
correlated, i.e. <f>cou(v 0 , • • • v n ) = 4>coii{v) n , then we have 
p( T 0i ■ ■ ■ T n \v 0 , ■■■v n ) = p(r ]u) n , yielding ( t 0 t u ) = ( r 2 ) 
and therefore C TT (n) = [(r 2 ) — ( r ) 2 ] 6 n , 0 - 

Now, let us make the hypothesis that free 
flight times are independent of the velocities, i.e. 
p(tq, ■ ■ ■ r n |i;o, • • • v„) = p(ro, ■ ■ ■ r n ). We want to discuss 
the implication of this simplification on the free path 
ACF. In general, since £ n = v n T n , the free path ACF is 
expressed as (£ 0 £n) = {T 0 T n v 0 v n ), where 


( T 0 T n v 0 v n ) = J dr 0 ■ ■ ■ dr n dv 0 ■ ■ ■ dv n T 0 T n v 0 v n 

p(TO, ■ ■ ■ T n \v 0 , ■ ■ ■ V n )(j) co ll(v 0, • • • V n ). 


( 21 ) 


Consistent to our assumption, we obtain (£ 0 £n) = 
(xoT n )(vov n ). Moreover, if we make the additional hy¬ 
pothesis that r n is an uncorrelated process, the free path 
ACF transforms to (£o£n) = i T ) 2 (v oV n ). Substituting 
this into Eq. (19) we then obtain: 


(dxl) 


(r ) 2 


n— 1 

n(v 2 ) +2 '^2(v 0 v m )(n-m) , 

m—1 


( 22 ) 


with (r) = /o°° drP{r )r and (v 2 ) = dv<fr coU (v)v 2 . 

The former approximation corresponds to a Poisson pro¬ 


cess, as we will discuss in Sec. [TV] We highlight that the 
only two assumptions made to define a Poisson process 
are i) independence between free flight times and veloc¬ 
ities in collisional representation, ii) no memory effects 
characterizing the free flight collisional process. The nu¬ 
merical evaluation of the approximated MSD in Eq. (22) 
is reported in Fig. H3a) and compared to the expression 
in Eq. (16 1 , showing a satisfactory agreement only in the 
low packing fraction limit. For more dense fluids, the in¬ 
termediate caging regime appears to be poorly captured 
by Eq. (22). The negative antipersistent part of the free 
path ACF (Fig. [9b, is indeed less pronounced in the ap¬ 
proximation (£ 0 Cn) — (t) 2 (v 0 V n ). 

We finally study the cross-correlation functions (£0 T n) 
and (ro£ n )- Both vanish since (£) = 0, we therefore con¬ 
sider the absolute values, i.e. (|£o|T n ) and (r 0 |£„|). Fig[l2] 
shows (|£ 0 |r„). The first term can be analyzed in the low 
packing fraction limit, reading: 

/ oo /»oo 

dv\v\(j) coU (v) / dTT 2 p(r\v) (23) 
-00 J 0 
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FIG. 11: (Color online) Poisson approximated expression for 
the MSD. (a) MSD in collisional representation for several 77 
calculated according to the Poisson approximation in Eq. (22 1 
(dashed lines). Th e so lid lines refer to the Green-Kubo ex¬ 
pression from Eq. (18 1 , coinciding with the Poissonian ap¬ 
proximation only in the low 77 limit, (b) MSD in continuous 
time representation using the first two terms of the Poisson 
formula appearing in Eq. <|39| ) (dashed lines), and the Green- 
Kubo expression in Eq. ( |15[ ) (solid lines). Here too the low 
77 systems appear to be well captured by the approximated 
formula. The high 77 limit requires a com plete determination 
of the correlations appearing in Eq. (38 1 . 


where we dropped the index 0 on the RHS, to simplify 
the notation. The conditional probability p(r\v) is expo¬ 
nential for low 77 (see Fig|2ja)), and substituting Eq. ([ 7 ]) 
we have: 


FIG. 12: Color online) Normalized cross correlation of free 
flight time and free path length. 




(l£o| r o> = 2(r) 2 r dv\v\^l. (24) 

J -00 Vcou(y) 

In ID (Jepsen gas), where tp C oii(v) = jjmbW, one ob- 
tains(fe|ro) = 2 3 / 2 (r) 2 y^j4l]. 


FIG. 13: (Color online) Intermediate scattering function, (a) 
Continuous time Fs(q,t) for several 77 calculated according 
to Eq.(|25[). The passage from the exponential decay to the 
stretched exponential regime, particular to caging in granular 
and colloidal systems is evident, (b) Collisional representation 
Fs(q,n), Eq. (26 1 . The exhibited trend traces that of panel 
(a). 
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D. Intermediate scattering function 


We now study a different dynamical function at the 
single particle level, the single component intermediate 
scattering function. In continuous time it is defined as: 

Fs(q,t) = (cos [q (x(t) — x(0))]). (25) 


This function is the real part of the x component of the 
single particle (incoherent) dynamic structure factor with 
a wave vector q, Fg(q,t) = (e lq, ( r ( t ) _r (°))) ; and has been 
the subject of an extremely large amount of numerical 
and theoretical work in recent years. While a single ex¬ 
ponential decay is expected for a normal fluid state, the 
emergence of a g-dependent 2 -step decay at high packing 
fractions has been recognized to characterize the behav¬ 
ior of systems in different areas of physics, ranging from 
supercooled liquids and granular materials to glasses and 
colloidal systems [701 ESI HUM]- Intuitively, this non¬ 
exponential regime (often referred to as a stretched expo¬ 
nential decay) should be related to the cage effect, which 
has been discussed as the main reason for the plateau-like 
region observed in the MSD (Fig.[7](a)) and the antiper- 
sistent back-scattering phenomenon exhibited by the ve¬ 
locity ACF (Fig. | 6 ](a)) . The 2 steps of the Fg relaxation 
dynamics consist in a fast and local “caged” /3-process, 
and a slow “cage escape” a-process, related to the cage 
restructuring, occurring in a more and more cooperative 
manner as the molecular packing increases [85]. The ap¬ 
parent stretching of the relaxation functions can also be 
explained by spatial heterogeneity, where relaxation oc¬ 
curs exponentially on different timescales in each spa¬ 
tial domain |j58J. In Fig. [l3](a) we report the numerical 
study of the intermediate scattering function Eq. (25) 
for several 77 . The progressive approach from the sin¬ 
gle exponential toward the 2 -steps stretched exponential 
relaxation is clearly shown, and, to our knowledge, this 
is the first time that such evidence is presented for a 
hard-disks system. This result, together with the plots 
in Figs. |t a ) ancl[F}a), demonstrates 2D hard-core sys¬ 
tems present several important analogies with complex 
colloids, glasses and supercooled liquids. 


Inspired by the previous analogy between continuous 
time and collisional dynamical quantities, we study the 
intermediate scattering function in collisional representa¬ 
tion, which takes the form: 


Fs(q, n) = (cos [q ( x n - x 0 )]). (26) 

Once again, if one neglects the early stage regime char¬ 
acterizing the ballistic behavior in the continuous time 
analysis, the decay exhibited by Fg(q,n ) in Figfl3](b) 
closely resembles that of Fg(q, t). In particular the expo¬ 
nential relaxation expected for the fluid state is recovered 
for low 77 , while the stretched exponential decay appears 
to be more marked as as 77 increases. 




FIG. 14: (Color online) Self part of the van Hove func¬ 
tion. (a) Continuous time Gs(x,t) for 77 = 0.3, calcu¬ 
lated according to Eq.( |27| ), for two different times t (solid 
lines). The dashed lines represent the Gaussian form 
,1_. = exp (— x 2 /2(5x 2 (t))), where the values of (5x 2 (t)) 

^2ir(Sx 2 (t)) 

have been drawn from FigQ a), (b) Collisional representa¬ 
tion Gs{x , n) for 77 = 0.3, calculated according to Eq.(28 1 , for 
two different values of n (solid lines). The dashed lines rep¬ 
resent the Gaussian form — 7 === exp (— x 2 /2{8x\)), where 


■\J2n(Sx '^ l ) 


the values of (<5a;„) have been drawn from Fig. p[a). 


E. Self part of the van Hove function 

The last dynamical quantity that we study at the single 
particle level is the self part of the van Hove function j99j. 
This function is the equilibrium, canonically averaged, 
probability distribution of the particle displacements (the 
propagator), that for technical reasons is used in experi¬ 
ments in its one-dimensional form: 

G s (x,t) = (6(x-[x(t.)-x{ 0)])). (27) 

Generally, in materials close to the glass or jamming 
transition and colloidal gelation, the non-Fickian char¬ 
acter of the single particle displacements is indicated 
by three factors: (i) the relaxation functions decay 
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FIG. 15: (Color online) Self part of the van Hove func¬ 
tion. (a) Continuous time Gs(x,t) for 77 = 0.695, cal¬ 
culated according to Eq. (271, for two different times t 
(solid lines). The dashed lines represent the Gaussian form 
= exp (— x 2 /2{Sx 2 (t))), where the values of ( Sx 2 (t )) 


y/2lT (Sx 2 (t)) 

have been drawn from Fig[7])a). Small but persistent devia¬ 
tions from the Gaussian expression in the tails of the distribu¬ 
tions are displayed at both times, (b) Collisional representa¬ 
tion Gs{x,n) for 77 = 0.3, calculated according to Eq.(28), for 
two different values of n (solid lines). The dashed lines rep- 

exp (—x 2 /2(5x 2 ) l )), where 


resent the Gaussian form 


yj2n(5x ^ l ) 


the values of (5x„) have been drawn from Fig[ 8 ])a). Deviations 
similar to the ontinuous time case appear. 


non-exponentially, (ii) the MSD exhibits a subdiffusive 
plateau at intermediate time scales, (in) the intermedi¬ 
ate scattering function Fg exhibits stretched-exponential 
decay, (iv) the self-part of the van Hove distribution func¬ 
tion is non-Gaussian. In the previous sections we have 
provided the numerical evidence of the first three phe¬ 
nomena in both continuous time and collisional repre¬ 
sentation, and we now tackle the last one. Moreover, 
the deviation from the Gaussian form is generally inter¬ 
preted as the key signature of the so-called dynamical 
heterogeneity 151] - he. when slow particles tend to 
cluster together forming cages, as well as faster parti¬ 


cles forming fluid-like regions [55U5CT1 f5Sl 1100)4102] . In 
Fig. |7~i|a) we report the van Hove function from Eq. (27) 


calculated for 77 = 0.56, in the dilute regime (solid 
line). The dashed lines represent the Gaussian expres¬ 


sion , 1 =e 2 <^ T U)) where the value of the MSD 

y/2n(5xHt)) 

has been obtained from data in Fig. [Tja). It is apparent 
that in this case the Fickian behavior is respected and the 
Gaussian form reproduces with remarkable accuracy the 
numerics at two different times. In Fig. |15[ a) the van 
Hove distribution function is displayed for 77 = 0.695, 
within the coexistence phase. Gs(x,t) has been calcu¬ 
lated at two different times according to the different 
stages attained by the MSD in Fig. [7](a): within the sub- 
diffusive plateau and the subsequent diffusive phase, once 
the tracer has escaped from the cage. In this case one 
can clearly see that the dashed lines accounting for the 
Gaussian Fickian behavior do not capture the tails of the 
van Hove function for times corresponding to the plateau 
region in Fig. [7](a). When the dynamics is diffusive, the 
Gaussian behavior is instead restored. The appearance of 
non-Gaussian tails is even more apparent for 77 = 0.713 
(Fig. [l 6 |(a)) where at short and larger times the expo¬ 
nential tails characterize the behavior of the van Hove 
function, suggesting that the particle dynamics can be 
represented by a number subsequent caging events. This 
is in line with the observation of plateau-like regions de¬ 
tected in the MSD (Fig.[7|a)), and with the velocity ACF 
antipersistent tails observed in Fig. [ 6 j)a). 

Studying the single component van Hove function in 
collisional representation through the definition 


G s (x, n) = (5 (x - [x n - x 0 ])) 


(28) 


traces the behavior exhibited in continuous time. In 
Figs. [ 3 b), [l^b) and dib) the solid lines represent 
the numerical data, whereas the dashed curves account 


for the Gaussian form 


= e 2 < 5 x D, where (<5x„) is 


^/‘2-k (dx 2 L ) 

drawn from the results in Fig. [9])a). For the lowest pack¬ 
ing fraction the agreement observed in continuous time 
is respected. For 77=0.695 shown in Fig.[l5|(b), the agree¬ 
ment between numerical curves and Gaussian expressions 
seems to be better for n = 1000 than n = 30. For 
77 = 0.713 (Fig.[l 6 |(b)) the exponential non-Gaussian tails 
characterize the subdiffusive regime as well as the ensuing 
diffusive phase. We also note the continuous-time prop¬ 
agator exhibits oscillations related to the disk packing, 
which seem to disappear in the collisional representation. 
This can be attributed to the fact that solid packing is a 
collective process, orchestrated with external time t. 


IV. FROM COLLISIONAL TO CONTINUOUS 
TIME REPRESENTATION 


We now wish to study the relation between continuous 
time and collisional representation. In Sec III C we have 
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provided evidence that the free flight time r n evolves ac¬ 
cording to a non-Markovian dynamics, and this property 
is fulfilled by gases at any 77 . The distribution of times 
t required for a particle to undergo n collisions, i.e. the 
probability that a particle makes exactly n collisions up 
to a time t, is given by: 

P(n, t) = J dr 0 - ■■ dr n -i dv 0 ■ ■ ■ dv n -1 S ^ r^j 

p(x Oj ' ' ' T~n— 11^0> * * * 'Un—l)4*coll (VO: * * * 'Cn— l)- 

(29) 

Hence, the exact expression of P(n,t) can be derived 
uniquely upon knowing the joint probability distribution 
of the velocities, and the ensuing conditional probability 
of the free flight times. Following Ref. M , we take the 
Laplace transform of both sides of Eq.(29): 


-*E" 


P(n, s) = / dr 0 ■ ■ ■ dr n -i dv 0 ■ ■ ■ dv n -ie 


p(r 0 , * ‘ 1 T n — 1 |t7o , * V n — 1 )<ficoll (vo j " * Vn— 1 ) • 


(30) 


Assuming the velocity process to be uncorrelated, the 
former Laplace transform then reads: 


P(n, s) = 


dvp(s\v)(f> coU (v) 


(31) 


where P(s) = f dv p(s\v)(f> co ii(v) is the Laplace trans¬ 
form of the free flight time distribution function P(r) 
in Eq. Q. In the dilute limit, the exponential form of 
p(r\v) leads to: 


P(n, s) = 


dv 


r )s4>coii{v ) + 4>mb(v) 


(32) 


where we made use of Eq. 0 The same expression has 
been derived in Ref. m assuming the velocity collisional 
process to be Markovian, where the numerical Laplace 
inversion of Eq. (321, or analogously of P(s) n , has been 
shown to closely agree with P(n,t) extracted from sim¬ 
ulations for 77 in the fluid phase. This tells us that the 
uncorrelation assumption of the velocity collisional pro¬ 
cess hold in the dilute limit, if one considers P(n,t) or 
P(t) [IB]. In the large packing fraction limit, however, 
this assumption fails and the inverse Laplace transform 
of Eq. (321 does not capture the probability P(n,t). 

Let us consider the case in which the free flight 
time distribution takes the exponential form P(t) = 
Aye"*/< r >. In view of Eq. Q, we are implicitly assuming 

p(r\v) = -^xe~ t P T \ i.e. a unique time scale character¬ 
izes the collisional process irrespective of the undergone 
velocity. This corresponds, de facto , to taking free flight 
times independent of velocities. We can thus write the 
former Eq. 0 as: 

P{n, t) = J dr 0 ■ ■ ■ dT n -x 6 ^t - ^2 t-^J P(r 0 , • • • r„_i), 

(33) 


and, if we make the further assumption of considering the 
free flight times to be uncorrelated, one finally obtains for 
the Laplace transform (301 : 


P{n,s) = 


1 


((t}s + 1Y 


(34) 


Inverting back in time domain, yields the Poisson formula 

G! : 


i / j. \ n—1 —ji-x 

1 / t \ e ” 


PM = V) \V) J (»—!)!■ 


(35) 


which has shown a large discrepancy to Eq. (29) even in 
the dilute limit [T31 Hi] . However, the collisional MSD 
expression in Eq. (22 1 derived in this framework seems 


to work fairly well in the low p limit (see the black curve 
corresponding to 77 = 0.3 in Fig{TT](a)) . 

We now turn to the average time between n collisions, 
i.e. the first moment of P(?z, f), within the context of 
the Poisson approximation, i.e. using the Poisson ap¬ 
proximated expression in Eq. ( |35[ ). A straightforward 
calculation yields: 


(t n ) = n(r ), 


(36) 


and, surprisingly, it reproduces very well the numerical 
values reported in Fig |17| not only for low 77 but also for 
77 in the solid regime. 

To conclude this section, let us turn to the single par¬ 
ticle propagator. To obtain the self part of the van 
Hove function in continuous time of Eq. (271 from 


the corresponding function in collisional representation 
of Eq. (28), we have to recur to the well-known subordi¬ 
nation relation m : 


G(x, t) = ^2 G(x, n)n(?r, t), 


(37) 


n—1 


where n(n, t) is the single particle probability of having n 
collisions at time t. The connection between the formerly 
introduced probability P(n,t) in Eq. (29) and n(n, t) is 
given, for large n, by n(n,f) = (r)P(n, t). Thus, the 
MSD can be easily cast ed in the following form: 


(Sx 2 (t)) = (r) ^ 2( 5x n) p ( n ’t )■ 


(38) 


n—1 


This expression constitutes the relation between the 
MSD in collisional representation and that in continu¬ 
ous time, but the rigorous evaluation at any 77 requires 
the knowledge of the conditional probability distributions 
in Eq. (29). However, a semianalytical expression of 


the MSD in continuous time can be achieved if we use 


the Poissonian approximation. Thanks to Eq. (22) and 


Eq. (35) we have, after elementary algebraic passages: 


{5x\t)) = (TY(v 2 ) + (r)(v 2 )t+ 


2EZ (v 0 v m )(n- m) f 

nr, O m — 1 \ \ / 


n—1 --At 

e vo 


(39) 


(n — 1)! 
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This approximated expression can be considered valid 
only in the low 77 regime and, as such, the three terms 
in it can easily be interpreted. The first term accounts 
for the ballistic regime, the second is the leading diffu¬ 
sive term which dominates the single particle dynamics 
after the first collision. The third term represents the 
logarithmic corrections expected for an hard-disks gas. 
It is interesting to notice that the second moment of the 
velocity ( v 2 ) is obtained using the collisional stationary 
distribution of velocity tp co ii(v), instead of the regular 
Maxwell-Boltzmann <Pmb(v). As a first approximation 
we neglect the logarithmic corrections in Eq. (39) which 
are believed to play a significant role for large times (see 
the discussion in Sec III B). In Fig. [TT|b ) we observe that 
the approximated Poisson formula ( |39[ ) accurately repro¬ 
duces the MSD in continuous time for low 77 , whereas it 
clearly fails to account for the high packing regime, as 
indeed is expected. 


V. DISCUSSION 

We have studied statistical and dynamical properties 
of a 2D system of hard-disks in collisional representa¬ 
tion, with packing fractions ranging from the fluid to the 
solid phase. Throughout this paper we tried to place our 
results withing a broader context, resulting in a compre¬ 
hensive study of 2D hard disk systems. The first part 
of our analysis focused on the statistics of velocity, free 
flight time and the path undergone by a disk between 
subsequent collisions. We have provided the numerical 
and theoretical evidence that the velocity x and y compo¬ 
nents are correlated in collisional representation, together 
with the fact that free flight times are strongly correlated 
to the on-collision velocities. This result has previously 
been investigated mmm only in the low packing frac¬ 
tion limit, and we have extended it to the coexistence and 
solid phase. We have shown that kinetic and Enskog the¬ 
ory do not furnish correct expressions for the average free 
flight time (r) and path (|£|) as a function of the packing 
fraction 77, not even in very dilute systems. Furthermore 
our numerical analysis spanning packing fractions up to 
the high 77 limit has unveiled plateau-like regions of (r) 
and (|£|) corresponding to the fluid-solid coexistence re¬ 
gion, which to the best of our knowledge have not been 
reported before. 

The second part of this study has been devoted to the 
dynamical single particle properties in collisional repre¬ 
sentation, and the comparison to their continuous time 
equivalents. We have shown a full analogy between ob¬ 
servables calculated in continuous time (as a function of 
t) and in collisional representation (as a function of the 
number of collisions, n), in spite of the fact that the pro¬ 
cess is not truly Poissonian [T3]. We have shown that the 
velocity is never a Markovian process both in continuous 
time and collisional representation. This is attributed to 
the persistent memory effects characterizing the velocity 
autocorrelation function at low packing fractions, that 


eventually turn into antipersistent tails at high packing 
fractions. In particular, this finding contradicts the com¬ 
mon assumption of Markovianity of velocity in the colli¬ 
sional representation at low packing fractions [15) . Fur¬ 
thermore we have analyzed the mean squared displace¬ 
ment, intermediate scattering function and self-part of 
the van Hove function, or propagator, both in continuous 
time and collisional representation, showing remarkable 
similarities with colloidal systems, glassy systems and su¬ 
percooled liquids close to the glass or jamming transition. 
These analogies provide a strong indication that within 
the coexistence phase the simplistic hard-disk model cap¬ 
tures properties of glassy behavior. 

Finally our analysis shows that, according to the ob¬ 
servable of interest, different approximations yield bet¬ 
ter descriptions of the data. Indeed we have shown 
that, although the system is generally non-Poissonian, 
at low packing fractions it can be fairly well described 
as a Poissonian model, i.e. considering uncorrelated free 
flight times independent of velocities. This approxima¬ 
tion works well if one considers the mean squared dis¬ 
placements ( S 2 x(t )) and ( S 2 x n ) or the average time as a 
function of the number of collisions ( t n ). On the other 
hand we have pointed out that a different approximation, 
considering the on-collision velocity process uncorrelated 
or Markovian but still correlated to the free flight times, 
provides excellent predictions in describing the proba¬ 
bility P(t) at low 77, the probability n(?7,,f) of having 
n collisions up to a time t in the dilute limit, or when 
one attempts to derive the on-collision stationary veloc¬ 
ity distribution function ipcoii(v) [13). However, none of 
the former assumptions is entirely correct, since velocities 
are correlated in collisional representation at any 77, and 
they are also correlated with free flight times: therefore 
none of them can furnish an overall satisfactory theoret¬ 
ical framework, not even in the very dilute regime. 
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Appendix A: Single component mean free-flight 
times relation 

Let us assume the sequence of velocities a particle ac¬ 
quires during N collisions is { 74 } where i £ [1, TV]. The 
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number of occurrences of some velocity v is denoted: 

N poo 


n , 


= £ 


dv8{ Vi — v ), 


(Al) 


and the total time spent on the velocity v is then given 
by T v = Y2? T i I dv6(vi — v). The velocity distribution 
function in collisional representation follows: 


n v 


Vcoiiiv) = 


(A2) 


and the velocity distribution function in continuous time 
can be approximated as: 


T 

( \ ± v 
<PMb(v) = WW. 

\ T ) 


(A3) 


Mul tipl ying and dividing the RHS by N, substituting 
Eq.(A2) and realizing that Ty can be approximated as 


T v = 71 v (t(v)) yields the relation: 

<Pmb(v) = ^ j' ^ (Pcoli(v). 

(v 


(A4) 


Appendix B: Mean free-path and mean free flight 
time from kinetic and Enskog theory 

The familiar expression for the average free-flight time 
arising from kinetic theory for a system of elastic hard 
disks is HD] : 


( T )kt = 


1 


4a/7 rksTno' 


(Bl) 


where n is the density number, i.e. n = dY. Recall¬ 
ing that the packing fraction is rj = nno 2 / 4, we obtain 

(r)fet = if Tj^k^r- Fig- |^J( a ) shows tliat E T ( B1 > agrees 

with the numerical data only qualitatively (blue line). 
The mean free-path is defined as: 


(l^l)fci — (| v\)cont {x)kt — 


707 

2 9/V 


(B2) 


The corresponding expression for the single component 
reads: 


(0 kt — (|^|) cont{x) kt — 


2? ! 2 r\ 


(B3) 


The mean free path value furnished by the Enskog theory 
for a system of 2D hard disks is given by: 


(l€l)o = 


2 9 / 2 rjY[r]\ : 


(B4) 


1/2 to match the kinetic theory expression in Eq. (B2|, 
as should be the case. The single component is then: 


(£)o = 


2 1/ 2 r]Y[r]) 


(B5) 


The expressions in Eqs. ( |B5| and (|B3j) are displayed in 
the inset of Fig. |4](b), showing a non-satisfactory agree¬ 
ment with the average free path obtained in the numer¬ 
ical simulations (red and blue lines respectively). The 
mean free flight time can be obtained through: 


\ r /o = 


(jo = g 

|v|) cont 16rjY[ri] \J k B T 


(B6) 


and is plotted in the inset of Fig. |4](a) with the expression 
in Eq. (B5) for the mean free path (red line). 


Appendix C: Generalized Kubo theorem in 
collisional representation 

Let us consider a discrete stochastic process such 
that (£„) = 0. The relation between the MSD and the 
symmetrized ACF is provided by Eq. (191, and brought 
again here: 


(&*£) = 2 £ C' K (m)(n - m). 


(Cl) 


m —0 

The on-collision transport coefficient D n is defined as: 

n—1 


°n=J2 Clid r 


(C2) 


m =0 


and it is possible to study its behavior once the depen¬ 
dence of C ££ on n is defined. Two main situations can 
arise: (i) D n ~ const , and ( ii ) D n ~ n 1- ' 3 (/3 > 0). 

(i) Suppose that is such that D n = 

V" n Ctt (m) ~ D with D = const. Hence from 
Eq (jcl} we have: 


(5xl) = 2D n - 2 £ 


(C3) 


m— 0 


The second term on the RHS of Eq. (C3) is 


Em=o C idm)m < Dn , yielding: 

{5x 2 n ) ~ 2D n 


(C4) 


for a large enough n. 

( ii ) Consider now a process for which C^(n) ~ c^n - ' 3 . 
The corresponding discrete transport coefficient is then 
given by 1104 :, 1051 : 


where Y[x\ = (l — y^x) /(I — x) 2 is the En skog factor. _ 

Notice that the expression furnished in Eq. (B4) differs 
from the formula provided in Refs. Earns] by a factor 


m 


(n-l) 1 - 

1-/3 


P 


dx- 


In— 1 


rP+i 


(C5) 
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where £(/3) is the Riemann zeta function and [•] repre¬ 
sents the integer part. Hence for large n we can write 

1-/3 

D n ~ D + cp n 1 _ /3 where D is a constant. From Eq. 
it then follows: 


Cl 


( Sx ~ 2 Dn + 2cp- 


(C6) 


(l-/3)(2 -13) 

If 0 < j3 < 1 the process is superdiffusive. Indeed cp is 


positive meaning that persistent memory effects charac¬ 
terize the free pat h Moreover D n diverges as ~ n 1 ^' 3 
according to Eq. (|C5|. If 1 < /3 < 2 and D = 0 the 


process is subdiffusive. In this case cp < 0, which defines 
the antipersistent tails of the symmetrized ACF C^(n). 
D n is positive and decays to zero as n 1- ^. If D ^ 0 the 
process is diffusive and D n ~ D. If/3>2,Z?^0, one 
has normal diffusion and D„ ~ D. 
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FIG. 16: (Color online) Self part of the van Hove func¬ 
tion. (a) Continuous time Gs(x,t) for r/ = 0.713, cal¬ 
culated according to Eq. (2TI, for two different times t 
(solid lines). The dashed lines represent the Gaussian form 
= exp (— x 2 /2{Sx 2 (t)}), where the values of ( Sx 2 (t )) 


y/2lT (Sx 2 (t)) 

have been drawn from Figj7|a). Deviations from the Gaus¬ 
sian form are clearly displayed: in particular, for the longer 
time corresponding to the asymptotic linear regime in Fig[7](a) 
non-Gaussian tails are reminiscent of what is found in granu¬ 
lar gases and colloids, (b) Collisional repr esen tation Gs{x,n) 
for r/ = 0.3, calculated according to Eq. (281, for two differ¬ 
ent values of n (solid lines). The dashed lines represent the 
Gaussian form . 1 . = exp (— x 2 /2{Sx%)), where the values 

W2tt(5x^ l ) 

of (Sxn) have been drawn from Fig. |8])a). Deviations from 
the Gaussian diffusive form appear. Moreover non-Gaussian 
tails similar to those in panel (a) characterize the free path 
dynamics, for n well beyond the plateau in Fig. [8](a) . We note 
that oscillations present in the continuous, representing spa¬ 
tial order due to high packing fractions, are not apparent in 
the collisional representation, due to the decoupling of space 
and time (detailed in main text). 
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FIG. 17: (Color online) Average time after n collisions, {t n ), 
for several r/ r/ (solid lines). Das hed lines correspond to the 
Poissonian approximation in Eq. (361. 
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FIG. 18: (Color online). Free path length probability dis¬ 
tribution (appearing in Fig. |3[b )) rescaled by: (a) the ki¬ 
netic theory expression (Eq. ( |B3| )) , (b) the average free path 
(FigMjb)) and (c) the free path furnished by the Enskog the¬ 
ory (Eq. None of them shows good collapse of the 

curves. 














